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: Avian coronavirus infectious bronchitis virus (IBV) poses economic threat to the poultry industry 

: worldwide. Pathogenic IBV 3575/08 was isolated from broilers vaccinated with the attenuated viral 

> vaccine derived from a Taiwan strain 2575/98. In this study, extensive investigations were conducted on 

: the genome sequences, antigenicity, pathogenicity, and host immune responses of several IBV strains 

: in specific-pathogen-free chickens. Sequence analyses revealed that 3575/08 and 2575/98 shared high 
homology in their structural genes, but not in non-structural accessory proteins such as 3a, 3b and 5b. 
Despite a high degree of homology in their spike protein genes, cross neutralization test showed low 

: cross protection between 3575/08 and 2575/98, suggesting distinct antigenicity for the two strains. 

: Animal challenge experiments exhibited strong respiratory and renal pathogenicity for 3575/08. In 

: addition, early and prolonged viral shedding and rapid viral dissemination were observed. Immune 

: gene expression profiling by PCR array showed chickens infected with 3575/08 had delayed expression 

: ofa subset of early innate immune genes, whereas chickens infected with the wild-type or attenuated- 

: type 2575/08 revealed quick gene induction and efficient virus control. In summary, this study reveals 

> anew IBV strain, which harbors a known local genotype but displays remarkably altered antigenicity, 
pathogenicity and host defenses. 


: Avian coronavirus infectious bronchitis virus (IBV) is a member of Coronaviridae family, and the virus is respon- 
- sible for infectious bronchitis, an acute, highly contagious chicken viral disease. Respiratory signs such as cough- 
: ing, sneezing, and nasal discharge are the primary symptoms associated with IBV infections, and some of the 
: strains also display kidney and oviduct tropisms, causing urogenital symptoms! *. The disease poses economic 
: threat to the poultry industry worldwide. More than 50 serotypes of IBV have been documented since the first 
: isolation of virus in Massachusetts (Mass), US in the 1930 s*. The most commonly detected genotype of IBV in the 
: US is Ark-DPI, accounting for 23% to 65% of total isolation per year*. While most strains circulating in the US 
: cause respiratory symptoms, nephropathogenic QX-like strains are dominant in Europe’. In the Far East, nephro- 
: pathogenic LX4, LDL as well as proventricular-type QX-like strains have been reported frequently’. Presently, 
: vaccination is adopted to manage the disease, and Mass type vaccine H120 is the most widely used vaccine for- 
: mulation”. Due to frequent strain variations, however, traditional Mass type vaccines can be ineffective against 
: infections by different IBV variants”>*. Effective disease management thus demands better elucidation of strain 
: variations. 

: In Taiwan, viral surveillance from 2005 to 2006 in poultry slaughterhouses showed the IBV flock prevalence 
: was 17% and the chicken prevalence ranged from 1.8 to 3.7%’. Most of the IBV isolates in Taiwan are nephrop- 
: athogenic strains*’. Also, IBVs in Taiwan are distinct from other genotypes in the world and can be divided into 
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two subgroups, Taiwan Group I (TW-I) and Taiwan Group II (TW-II)'*!!. As TW-I strains begin to emerge in 
neighboring areas in recent years!” "', there is a rising urgency to thoroughly investigate the TW IBV strains. 

IBV virion has a linear, positive sense and single-stranded genomic RNA. Its 27.6kb genome contains at least 
ten open reading frames (ORFs), 5’- 1a-lab-S-3a-3b-E-M-5a-5b-N-3’, encoding four structural proteins and sev- 
eral non-structural proteins (nsp)*. The four structural proteins include spike glycoprotein (S), envelope protein 
(E), membrane glycoprotein (M), and nucleocapsid protein (N). Of the different structural proteins, the S protein, 
which comprises of $1 and S2 subunits, is responsible for host cell attachment and virus-cell membrane fusion 
and is the primary antigen target for virus neutralization. Currently, genetic grouping for IBV is based primarily 
on the S1 gene’*. Among the genes that encode for non-structural proteins, the ORF 1a and lab encode two 
polyproteins, which are further cleaved into 15 nsps. These and other non-structural accessory proteins 3a, 3b, 
5a and 5b are associated with RNA replication, virulence, and host-pathogen interactions of IBV>. However, the 
functions of these nsps remain largely unknown and warrant further studies. 

Host innate antiviral immunity plays a critical role in viral clearance. In general, toll like receptor (TLR) fam- 
ilies, retinoic acid-inducible gene I (RIG-I) and melanoma differentiation-associated gene 5 (MDAS5, encoded by 
IFIH1 gene) are intracellular pattern-recognition receptors (PRRs) that function as viral infection detectors and 
all contribute to type I interferon (IFN) production’’. As RIG-I is absent in chickens, IFIH1 is a functional com- 
pensate'®, Avian TLRs are different from their mammalian counterparts, lacking TLR9 but possessing TLRILA, 
poultry-specific TLR15, and a TLR9 homologue TLR21’”. TLR3 and TLR7 are known to be involved in RNA 
virus recognition, and its stimulation leads to activation of IRF3 and IRF7 (IFN-regulatory factor 3 and 7). These 
factors then signal IFN transcription to trigger IFN-stimulated genes (ISGs) through JAK-STAT pathway to 
express various antiviral proteins'>!”'*. A global gene expression profiling of chicken kidney has revealed that 
103 immune and inflammatory-related genes are involved in the host defense response during the first week of 
IBV infection’. The TLR3/7 and MDAS signaling pathways and innate immune cytokine were induced after IBV 
infection in chicken tracheas at 3-8 days post IBV M41 infection”. In addition, STAT1 expression has also been 
found to be significantly elevated in response to IBV infection in ovo”!. However, there is a lack of understanding 
of early antiviral immune responses induced by various IBV strains, particularly responses triggered by strains 
with close genetic relationship but distinct pathogenicity. 

To control IBV in Taiwan, an embryo-attenuated strain (2575 AT) derived from IBV TW-I strain 2575/98 
wild type (2575 WT) has been developed as a vaccine®. These two strains were previously investigated to have 
close antigenicity* and limited sequence changes”. However, a novel TW-I strain 3575/08 emerged from a broiler 
chicken farm where the 2575 AT vaccine was applied. IBV 3575/08 has been observed to cause high morbidity 
and mortality in affected broilers. In the present study, we compared the viral genome sequences of IBV 3575/08 
and 2575/98 and investigated the infectivity and host immunity induced by 3575/08, 2575 WTT, and 2575 AT 
infections. To this end, sequence analysis, neutralization test, experimental infection, and immune-related gene 
profiling were performed. The results show that IBV 3575/08, which possesses heterologous ORF 3a, 3b, and 5b 
sequences to 2575/98, caused a more deadly infection with early and prolonged viral shedding, stronger clinical 
and pathological manifestation, and delayed innate immune responses as compared to 2575 WT and 2575 AT. 


Results 

3575/08 and 2575 WT share sequence homology except in 3a, 3b, and 5b genes. The complete 
genome of IBV 3575/08 was sequenced in this study, and the sequences were submitted to GenBank with the 
accession number KX266757. Excluding the 5’ and 3’ untranslated regions the entire genome size is 26,704 
nucleotides. The genome comprises ten ORFs in the order of 5/-1a-1lab-S-3a-3b-E-M-5a-5b-N-3’. No significant 
protein truncation was observed. To evaluate the nucleotide sequence identities among 3575/08 and other refer- 
ence IBV strains, sequence alignment and phylogenetic analyses on each gene segment were performed (Fig. 1 
and Supplementary Fig. $1). Comparisons of the available full-length genome showed that the sequence identi- 
ties between 3575/08 and other IBV strains ranged from 85.0% (with strain BJ) to 94.9% (with strain 2575/98). 
Overall, IBV 3575/08 showed close genetic grouping with 2575/98 and other Taiwan local strains in 1a, 1b, S, E, 
M, 5a and N genes. The sequence similarities of these protein genes between 3575/08 and 2575/98 were 95.6% 
(1a), 94.2% (1b), 95.6% (S), 93.3% (E), 92.4% (M), 93.4% (5a), and 96.3% (N) for the respective genes. Further 
comparison of the two spike protein subunits between 3575/08 and 2575/98 showed sequence similarities of 
97.1% and 94.4% in the S1 and S2 subunits, respectively. Based on the S gene, 3575/08 was clearly defined as a 
TW-I strain. The sequence identities of S gene between 3575/08 and other analyzed strains ranged from 67.6% 
(with strain Georgia 1998 pass8) to 96.4% (with strain 3468/07). While 3575/08 was closely related to 2575/98 
in four structural protein genes, their homology in ORF 3a (86.1%), 3b (82.1%), and 5b (95.5%) were relatively 
low (4.5% difference is considered significant for the 5b protein since its variation among IBV strains is small). 
According to the phylogenic tree (Fig. 1), IBV 3575/08 and 2575 WT are grouped in different clusters with respect 
to their 3a, 3b and 5b genes. 


Identification of IBV 3575/08 as a unique serotype. To investigate the antigenic serotype of the newly 
emerged TW-I strain 3575/08, cross (two-way) neutralization tests against three major serotypes in Taiwan— 
TW-I, TW-II and Mass type— were performed in specific-pathogen-free (SPF) chicken embryos. The neu- 
tralizing titers and the r-values were then analyzed (Table 1 and Supplementary Table S1). With the r-value of 
homologous virus-serum set as 100, the r-value directly correlates with the relatedness of the strains. The results 
showed that all the r-values obtained in this study were below 30, indicating all the tested strains belonged to dif- 
ferent serotypes. Particularly, TW-I 3575/08 conferred poor cross neutralization not only to heterogeneous Mass 
type H120 (r-value= 1.1) or TW-II 2296/95 (r-value = 3.4) but also to homologous TW-I 2575/98 (r-value = 3.3), 
supporting IBV 3575/08 as a unique serotype. 
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Figure 1. Phylogenetic analyses of 1a, 1b, S, 3a, 3b, 5a and 5b protein genes. The phylogenetic trees were 
constructed using the MEGA software version 4 by the neighbor-joining method (bootstrapping for 1,000 
replicates with its value >70%). The 3575/08 and 2575/98 strains analyzed in this study are indicated with red 
triangles. (AU: Australia, CN: China, NL: Netherlands, TW: Taiwan, UK: United Kingdom, US: United States, 
WA: Western Africa) 


IBV 3575/08 showed enhanced pathogenicity as compared to 2575 WT and 2575AT in chickens. 
To investigate the pathogenicity of IBV 3575/08, 2575/98 WT and 2575 AT, SPF chickens were intranasally 
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IBV strain | 3575/08 (TW-1) | 2575/98 (TW-1) | 2296/95 (TW-II) | 120 (Mass) | 
3575/08 (TW-I) 100 | 

2575/98 (TW-1) 33 100 

2296/95 (TW-II) 34 DR 100 

11120 (Mass) 7 26 D8 100 


Table 1. R-value of cross neutralization test with TW-I strains (3575/08, 2575/98), TW-II strain (2296/95) 
and Mass type (H120). 


challenged with 10° EIDs» of the three strains and monitored for survival and clinical signs on a daily basis. In 
1-day-old chicks during a 21-day observation period, IBV 3575/08 infection resulted in the lowest survival rate 
(16.7%, 1/6) as compared to 2575 WT (66.7%, 4/6), 2575 AT (100%, 6/6), and PBS control (100%, 6/6) (Fig. 2A). 
The earliest death from 3575/08 or 2575 WT-infected group was found on day 4 post infection (dpi). 

In addition, clinicopathological assessment for the three IBV strains was performed. Following intranasal 
challenge in 2-week-old chickens with 10° EID.) of the different IBV strains, the clinical scores (Fig. 2B,C) were 
recorded and analyzed (n=5 per group, see Methods for detailed clinical score interpretation). It was observed 
that beginning at 2 dpi, chickens infected with 2575 WT exhibited moderate tracheal rales (score = 1), persisting 
through 3 dpi (score = 1) and resolved at 4dpi (score = 0.2) (Fig. 2B). IBV 3575/08-infected chickens showed 
similar but more severe clinical signs at 2 dpi (score = 1.4) and 3 dpi (score = 1.2). At 4dpi, more than half of 
3575/08-infected chickens still suffered from mild respiratory symptoms (score = 0.6), and clinical signs receded 
after 5 dpi. In contrast, IBV 2575 AT-infected and PBS control groups did not show any clinical manifestations 
throughout the observation period. Collectively, over a 7-day observation period, the averaged clinical scores of 
IBV 3575/08-, 2575 WT-, 2575 AT-infected chickens were 0.49, 0.32, and 0, respectively (Fig. 2C). Overall, IBV 
3575/08 infection caused stronger clinical manifestations in chickens than 2575 WT and 2575 AT. 

The pathological lesions of experimentally infected chickens were further evaluated at 7 dpi, and representa- 
tive pictures are shown in Fig. 3. In the trachea epithelium, normal goblet cells and ciliated cells were seen in the 
PBS control group (Fig. 3D). However, the IBV 3575/08 infected chickens presented pathological lesions such 
as desquamation of ciliated cells, infiltration of heterophils, and lymphocytes and epithelial hyperplasia in the 
trachea (Fig. 3A, arrowed). Lesions in the 2575 WT (Fig. 3B) and 2575 AT (Fig. 3C) groups were similar but less 
severe as compared to the 3575/08 group. Upon observing the kidneys, the control group showed normal renal 
tubular epithelial cell with homogeneous plasma and distinctive renal tubular cavities (Fig. 3H). In the kidneys 
of IBV 3575/08-infected chickens, however, renal tubule dilation, necrotic tubular cells, and urate in the tubular 
cavity were observed (Fig. 3E, arrowed). Vesicular degeneration was also observed in the renal tubular epithelial 
cells from the 2575 WT group (Fig. 3F, arrowed). On the other hand, the 2575 AT group showed normal kidney 
patterns similar to the control group (Fig. 3G). The pathological findings indicate that both IBV 3575/08 and 
2575 WT are respiratory and nephropathogenic strains. IBV 3575/08 showed higher virulence as it inflicted more 
severe trachea and renal tissue damages in chickens than the 2575 WT strain. No virulence in the kidney was 
observed for IBV 2575 AT. 


More rapid and prolonged viral dissemination of IBV 3575/08 in chickens. To further investigate 
the viral dissemination of IBV 3575/08, 2575 WT, and 2575 AT in chickens, swabs and tissues were harvested 
from virus-infected chickens (10° EID.o, intranasal) at various time points. Viral detection was carried out by 
reverse transcription polymerase chain reaction (RT-PCR). As indicated in Table 2, 100% of the throat swabs 
from all groups were detected positive at 1, 2, 3, 7, 14dpi. At 21 dpi, both the 2575 WT and 2575 AT groups 
showed a decreased positive rate at 66.7% (2/3), whereas a 100% (3/3) positive rate remained for the 3575/08 
group. More importantly, for the cloacal swabs, virus was detected from 3575/08-infected chickens as early as 
1 dpi (60%, 3/5), whereas the earliest time points for 2575 WT and 2575 AT were 2 dpi (40%, 2/5) and 3 dpi (20%, 
1/5), respectively. There was no detectable level of viruses in any of the swab samples from the PBS control. 

To examine viral distributions, the lungs and kidneys, two major organs associated with IBV tropism, were 
harvested at 3, 7, 21 dpi. As indicated in Table 3, whereas all of the infected chickens from the 3575/08 or 2575 
WT groups were detected positive in the lungs at 3 and 7 dpi, only 2/3 and 1/2 chickens were detected positive in 
the 2575 AT group. As for the kidneys, virus was detected from the 3575/08-infected chickens as early as 3 dpi, 
whereas no virus was detected from the 2575 WT group until 7 dpi. In addition, IBV 3575/08 was detected in the 
trachea, proventriculus, bursa, and oviduct from experimentally infected chickens at 7 dpi (data not shown). In 
the 2575 AT group, no animal was detected positive in the kidneys across the indicated time points. By 21 dpi, 
virus was cleared from the lungs and kidneys in every group. These results show that IBV 3575/08 causes more 
rapid and prolonged dissemination and shedding in chickens as compared to the other two strains. The results 
also highlight the nephropathogenic nature of 3575/08 and 2575 WT in chickens. 


3575/08 infection elicited a similar level of antibody titer as 2575 WT in chickens. To examine 
the antibody response elicited by IBV infection in chickens, sera were collected from the 3575/08-, 2575 WT-, and 
2575 AT-infected groups at 14 and 21 dpi. Antibody titers were analyzed by an IBV antigen coated-ELISA, and 
3575/08 and 2575 WT infections elicited similar mean titer levels of 6,547 and 6,870 at 21 dpi respectively (Fig. 4). 
The antibody titer elicited by 2575 AT infection was lower (3,294) at 21 dpi. Similar titer results were obtained 
from the different infection samples at 14 dpi. 
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Figure 2. Percentage survival and clinical score of IBV infected chickens. (A) Survival rate of 1-day-old 
chicks after intranasal infection with IBV 3575/08, 2575 WT, 2575 AT. Chickens receiving PBS were served 
as control. (n= 6 per group) (B) Following intranasal infections in 2-week-old chickens, the daily averages of 
observed clinical scores within 6 days were recorded. Data was presented as means +s.e.m. (n=5 per group) 
(C) Mean clinical scores of each group during the observation period are indicated (n=5 per group). See 
Methods for detailed clinical score interpretation. ns: non-significant. **P < 0.01. ***P < 0.001. 


Antiviral immune response was delayed in IBV 3575/08-infected chickens. Next, to dissect the 
antiviral gene induction in chickens in response to IBV 3575/08, 2575 WT and 2575 AT infections, we examined 
a panel of immune response gene expression in the lung at 24 and 48 hours after the IBV infections using a quan- 
titative PCR array kit. The fold changes in the expression of 84 genes from infected chickens were normalized to 
that of naive chickens (Fig. 5A, also see Supplementary Table S3 for statistical data and Supplementary Fig. $2 
for the Venn diagram). At 24 hours post infection (hpi), 20 (23.8%) and 33 (39.3%) of the genes were induced by 
more than 3-fold over the naive control in the 2575 WT- and 2575 AT-infected chickens respectively, whereas 
only 5 (6%) of the 84 genes were upregulated in the 3575/08-infected chickens. At 48 hpi, 16 (19%), 3 (3.6%), 
and 3 (3.6%) genes were induced by more than 3-fold over the naive group in the 3575/08-, 2575 WT-, and 2575 
AT-infected chicken, respectively. Overall, a robust immune-related gene expression was rapidly induced in the 
2575 WT- and 2575 AT-infected chickens at 24 hpi, and a 90% overlap was observed among the upregulated genes 
between the two groups. These results identified the critical components in early host defense against IBV 2575 
WT and 2575 AT. At the same time point, 3575/08-infected chickens showed heightened expression of only 5 
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Figure 3. Pathological examination of trachea and kidney from IBV infected chickens. The pathological 
lesions of experimentally infected chickens were evaluated at 7 dpi, and representative pictures were shown. 
Panels D and H refer to control tissues. IBV 3575/08 infected chickens presented pathological lesions such as 
desquamation of ciliated cells, infiltration of heterophils, and lymphocytes and epithelial hyperplasia in the 
trachea (A, arrowed). The lesions in 2575 WT (B) and 2575 AT (C) groups were similar to the 3575/08 group 
but milder. For IBV 3575/08-infected chicken kidney, renal tubule dilation, necrotic tubular cells, and urate in 
the tubular cavity were observed (E, arrowed). Vesicular degeneration was found in renal tubular epithelial cells 
from the 2575 WT group (F, arrowed). For the 2575 AT group, the renal cells were normal as the control group 
(G). 


innate immunity genes (TLR1LA, TLR3, IRF7, STAT1, and C3) (Fig. 5B). 6 other genes upregulated by 2575 WT 
and AT at 24hpi, including TLR4, MyD88, IFNA3, IFNGR1, CD28, and CATH2, showed delayed stimulation fol- 
lowing the 3575/08 infection (Fig. 5C). In addition, 10 other genes that were not enhanced in either the 2575 WT 
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3575/08 100% (5/5) 100% (5/5) 100% (5/5) 100% (5/5) 100% (3/3) 100% (3/3) 


2575 WT | 100% (5/5) | 100% (5/5) | 100% (5/5) | 100% (5/5) | 100% (3/3) | 66.7% (2/3) 
2575 AT 100% (5/5) | 100% (5/5) | 100% (5/5) | 100% (5/5) | 100% (3/3) | 66.7% (2/3) 
Cloacal swab % IBV positive 


3575/08 60% (3/5) 60% (3/5) | 100% (5/5) | 100% (5/5) | 100% (3/3) | 100% (3/3) 
2575 WT 0% (0/5) 40% (2/5) | 100% (5/5) | 100% (5/5) | 100% (3/3) | 100% (3/3) 
2575 AT 0% (0/5) 0% (0/5) 20% (1/5) 80% (4/5) | 100% (3/3) | 100% (3/3) 


Table 2. Rate of IBV detection by RT-PCR in chicken throat and cloacal swabs after infections. 


3575/08 3/3 2/2 0/3 1/3 2/2 0/3 
2575 WT 3/3 2/2 0/3 0/3 2/2 0/3 
2575 AT 2/3 1/2 0/3 0/3 0/2 0/3 


Table 3. IBV detection by RT-PCR in chicken lungs and kidneys at 3, 7, and 21 day post infection. 


titer 


3575/08 2575 WT 2575 AT 
IBV strain 


Figure 4. Serum antibodies of IBV infected chickens. Chicken serum was collected at 14 and 21 dpi, 
respectively. ELISA was performed to evaluate the antibody level against IBV. Bars represent means + s.e.m. 
(n=5 per group). 


or the 2575 AT group were found to be upregulated in the 3575/08 group at 48 hpi. The finding suggests that IBV 
3575/08 may have evaded the initial virus sensing and delayed the signaling pathways for early antiviral defense. 
In 2575 WT- infected chickens, the genes upregulated at 48hpi (TLRILA, TLR3, and CATH2) overlapped with 
those at 24hpi. Similarly, the genes upregulated by 2575 AT at 48 hpi (TLRILA, IFNA3, and CD28) also over- 
lapped with those observed at the earlier time point. In addition, quantification of TLR3 and TLR4 expression 
levels by regular real-time RT-PCR was carried out to validate the results from the PCR array. The expression 
patterns of these two genes were consistent between the real-time PCR and the PCR array analysis, showing acti- 
vation of TLR3 but not TLR4 level at 24 hpi in the 3575/08 group (Fig. 5D). 


Discussion 

The spread of TW type IBV in Taiwan has been attributed to poor protection by commercial Mass type vaccines”’. 
An effective TW type vaccine strain 2575 AT has been previously developed through serial passages in SPF 
embryonated eggs to control IBV infection in Taiwan. However, a new IBV strain 3575/08 emerged from 2575 
AT-vaccinated broiler chickens and caused high morbidity and mortality. To characterize the new strain 3575/08, 
comparative analyses on viral genotype, serotype and pathogenicity were conducted. 

Pairwise sequence comparison and phylogenetic analysis revealed that 3575/08 and 2575/98 share high 
homology in most of the structural protein genes, including the S1 gene that are commonly applied for IBV strain 
genotyping. Therefore, IBV 3575/08 was grouped in the TW-I genotype with 2575/98 and other strains. It has 
been reported that the serotype of IBV is associated with the hypervariable region (HVR) 1 in the N-terminus 
(residues 45-114) of the S1 gene”, which overlaps with one of the antigenic sites located in the spike protein of 
IBV”. Despite high sequence similarities between 3575/08 and 2575/98 were found in the spike protein gene 
(95.6%), S1 subunit gene (97.1%) and HVR 1 (94.8%), a distinct serotype of 3575/08 was clearly established in 
this study through a cross neutralization test. Similar findings were previously reported in UK showing that IBV 
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Figure 5. Comparison of immune related genes in IBV infected chickens. (A) Profiles of immune related 
gene expression in the lungs after infection by different strains as evaluated by RT? Profiler Chicken Innate & 
Adaptive Immune Responses PCR Array at 24 and 48 hpi. Gene expression of transcripts with a difference of 
more than 3-fold over naive was illustrated (black, upregulated; gray, no change; white, downregulated). 

(B) Gene induction in the lungs of IBV 3575/08, 2575 WTT, and 2575 AT-infected chickens at 24 hpi was 
expressed as fold-increase over naive. (C) Gene induction in the lungs of IBV 3575/08 and 2575 WT-infected 
chickens was expressed as fold-increase over naive at 24 or 48 hpi. (D) Expression level of TLR3 and TLR4 

at 24 hpi in was validated by real-time RT-PCR. Fold changes were compared by ANOVA followed by 
Dunnett multiple comparisons tests. Dotted lines denote a 3-fold change over naive chickens. Bars represent 
means +s.e.m. (n=3 per group). *P < 0.05, **P<0.01. 
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strains of three distinct serotypes possess as little as 2% of amino acid variations in region 19-122 and 251-347 
of the S1 subunit”®. The variations led to failure of cross-protection”’, supporting that small sequence differences 
can result in a serotype change. In addition, it has been reported that small sequence changes in the S2 subunit 
may alter interactions with the $1 subunit, which could in turn affect the conformation of the S1 subunit and 
serotype-specific epitopes”*. A previous study on murine hepatitis virus (MHV), a group II coronavirus, also 
showed that a single amino acid substitution in region 1109-1116 of S2 subunit can confer resistance to neu- 
tralization by $1 subunit-specific monoclonal antibody”. Therefore, the gene variations in the $2 subunit (5.6%) 
between 3575/08 and 2575/98 may have a significant contribution to the serotype difference. Future studies are 
warranted to further examine the determinants of antigenicity. 

In addition to the new serotype, we also found that IBV 3575/08 exhibited higher mortality and more 
profound clinical signs as compared to 2575 WT and 2575 AT through experimental challenges in chickens. 
Following intranasal infection, IBV 3575/08 had a very short latent period. The viral shedding from cloaca was 
detected as early as 24 hpi and lasted for 21 days, and a rapid viral dissemination to renal tissues (3 dpi) was 
observed. As it has been documented that spike proteins play a major role in determining viruses’ cell tropism 
and pathogenicity”, the increased IBV 3575/08 virulence is a notable finding given the high degree of homology 
between the spike protein genes between IBV 3575/08 and 2575 WT. 

A previous study has shown that different IBV strains can cause differing IFN levels at various organs and 
time points in chickens*!. We hypothesized that a different innate immune response may account for the altered 
pathogenicity of IBV 3575/08. Therefore, the early immune-related gene expression in the lung following 3575/08, 
2575 WT and 2575 AT infections in chickens was analyzed by PCR arrays. In contrast to the 20 and 33 genes 
activated in the 2575 WT- and 2575 AT-infected chickens, only a small subset of the genes was triggered in the 
3575/08-infected chickens at 24 hpi, and a delayed activation was observed at 48 hpi. This result suggests that IBV 
3575/08 has the potential to evade the initial virus sensing and suppress the first wave of host defense, leading 
to enhanced pathogenicity observed in this study. Many coronaviruses have the ability to evade detection by the 
host’s immune system and disrupt the host immune response through multiple mechanisms”. It has been 
reported that small accessary proteins 3a and 3b of IBV can delay the onset of IFN response by interfering with 
the PRR sensing without disrupting the signaling pathways*. Previous studies have also demonstrated that 3a 
and 3b accessory proteins of IBV are not essential for replication but can modulate the IFN-B response at the 
transcriptional and translational level***”. In addition, 5b accessory protein of IBV has been found to shut off 
host immune response and limit interferon production™. It has been shown that a natural recombinant strain 
containing partial nsp regions from a low pathogenic H120 strain in the frame of a nephropathogenic strain 
possesses attenuated pathogenicity, which suggests strong associations between nsps and IBV virulence”. Studies 
on other coronaviruses such as MHV and severe acute respiratory syndrome coronavirus (SARS-CoV) have also 
shown that nsps are associated with viral replication and evasion of host innate immune response via viral RNA 
modification’”“!, host mRNA degradation, and IFN antagonist competition”. Therefore, given the low sequence 
homology in the accessory proteins 3a, 3b, and 5b between IBV 3575/08 and 2575 WT, we conjecture that these 
accessary proteins contribute to the blockade of PRR expressions, delayed immune response, and enhanced 
virulence of IBV 3575/08 in the present study. Future studies are warranted to further dissect the roles of IBV 
3575/08’s nsps and accessary proteins in modulating host’s antiviral defenses. An in vitro infection model with 
detailed functional analysis of nsps would shed clearer insights on the host-virus interactions. 

Also of note is that the number of PRR genes activated at the early time point correlated with IBV virulence in 
the study. For instance, IBV 2575 AT infection activated multiple PRR genes, including IFIH1, TLRILA, TLR2-2, 
TLR3, TLR4, TLR7, and TLRI5. In addition to the two primary sensors for IBV, TLR3 and TLR7*, other PRR 
molecules have been studied for their roles in coronavirus infections. In previous studies, SARS-CoV and MHV 
infections have been shown to trigger IFIH1 expression for MDAS signaling**”. It has also been reported that 
TLR4 expression can be induced by SARS-CoV infection“. The robust expression of PRR molecules triggered by 
IBV 2575 AT sheds light on the molecular mechanism behind the strain’s attenuated virulence and pathogenicity. 
In contrast, the lower number of stimulated PRR genes triggered by IBV 3575/08 and 2575 WT serve to explain 
these viruses’ increased virulence. 

The present study may also offer insight on IBV disease management and vaccine design. Firstly, identifi- 
cation of a new pathogenic IBV strain with a unique serotype and strong host tropism suggests a rapid viral 
mutation rate in chicken hosts. This observation raises the need for regular vaccine update to keep up with the 
mutating viral strains. As the antigenicity or pathogenicity shift may result from genetic changes outside of struc- 
tural protein genes, antigen derived from a single epitope (eg. spike protein) may not be sufficient in eliciting 
cross-protection against different serotypes. Multi-epitope vaccine formulations should be considered for better 
disease management. Secondly, the present study elucidates several key innate immune genes and TLR pathways 
(eg. TLR3, TLR4, MYD88, IRF7, IFNA3, etc.) responsible for early restriction of IBV infection. When designing 
an inactivated or subunit vaccine, immunomodulatory adjuvants may be tailored to promote an infection-like 
immune stimulation as observed in the study. Triggering innate immune gene expressions and TLR pathways in 
a virus-like manner may facilitate a robust adaptive immune response against IBV infection. 

In summary, this study compares the sequence, antigenicity, pathogenicity, and immune host defenses among 
IBV strains. The strain 3575/08, which possesses heterologous ORF 3a, 3b, and 5b sequences, exhibits a distinct 
serotype and strong clinicopathological manifestation. This virus can also evade the initial virus sensing and delay 
several early innate immune genes, collectively leading to enhanced pathogenicity in a chicken host. 


Methods 

Ethics statement. Experiments involving chickens were conducted at National Taiwan University, under 
an approved Institutional Animal Care and Use Committee (IACUC) protocol (no. NTU-103-EL-3). All animal 
experiments were carried out in accordance with the approved guidelines. 
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Virus. IBV strain 2575/98 was isolated in Changhua, Taiwan in 1998°, and designated 2575 wild-type (2575 
WT) throughout this study. Embryo passage 74 of 2575/98 was obtained as previously described’, and designated 
2575 attenuated-type (2575 AT) throughout this study. IBV strain 2296/95 was isolated in Taoyuan, Taiwan in 
19958. Mass type strain H120 was purchased from Abic Biological Lab (Beit Shemesh, Israel). IBV strain 3575/08 
was isolated from a broiler farm in Yunlin, Taiwan in October 2008. In that affected farm, death was first observed 
in 25-day-old broilers with typical IB disease symptoms. The morbidity was higher than 15%, and the mortal- 
ity was approximately 7.5%. All IBVs used in this study were propagated in 10-day-old SPF embryonated eggs 
(Animal Health Research Institute, Tamsui, Taiwan) via an allantoic route as previously described’. All the used 
viruses were ensured free of other avian respiratory pathogens, ie. avian influenza virus and Newcastle diseases 
virus contamination by RT-PCR“. IBV viral titer (egg infective doses 50%, EIDs9) was evaluated in 10-day-old 
SPF eggs and calculated by Reed-Muench method”. 


Viral genome sequencing and sequence analysis. Viral RNA was extracted from virus-infected allan- 
toic fluid using a viral nucleic acid extraction kit (Geneaid Biotech Ltd., Taipei, Taiwan) following the manufactur- 
er’s protocol. Reverse transcription was carried out with M-MLV reverse transcriptase (Invitrogen, Carlsbad, CA), 
and genome segments were amplified by PCR using primers listed in Supplementary Table $3. The PCR reaction 
mixture contains 0.5 pl of Phire hot start I] DNA polymerase (Thermo, Carlsbad, CA), 5 yl of 5X buffer, 0.5 pl 
of 501M forward primer and reverse primer (Tri-I Biotech, Taipei, Taiwan), 3 pl of 2.5mM dNTPs (GeneTeks 
BioScience, Taipei, Taiwan), 5 ul of CDNA and 10.5 il of nuclease-free water. For PCR, the initial denaturation was 
at 95°C for 3 min, followed by 35 cycles of denaturation at 95°C for 1 min, annealing at 55°C for 30s, and exten- 
sion at 72°C for 2 min. The final extension was at 72°C for 10 min. All the amplified PCR products were analyzed 
in 1% agarose gels electrophoresis. DNA sequencing was conducted by a commercial service (Tri-I Biotech). Each 
nucleotide was determined from at least three identical results generated from separate PCR products. Nucleotide 
sequences were assembled using DNAStar software (DNAStar, Madison, WI). Sequences of the reference IBV 
strains were retrieved from GenBank with the accession numbers listed in Supplementary Table S4. Multiple 
sequence alignments and sequence identity were analyzed with the Clustal W method available in BioEdit 7.2.5 
software. Phylogenetic trees were constructed with the neighbor-joining method using MEGA software version 
4°°, and the bootstrap values were determined from 1,000 replicates of the original data. 


Antisera production and cross neutralization test. Antisera against 3575/08, 2575/98 and 2296/95 
were prepared in 3-week-old SPF chickens as previously described™. The antiserum against Mass-type H120 was 
purchased from Charles River Laboratories (North Franklin, CT). Cross neutralization test was performed with 
the protocol previously described**. Briefly, four-fold diluted sera were incubated with the same volume of 100 
EID; homologous or heterologous virus at room temperature for an hour, respectively. The virus-serum mixtures 
were then inoculated into the allantoic cavity of 10-day-old SPF embryonated eggs. Seven days after inoculation, 
the eggs were opened and examined for typical lesions caused by IBV infection (embryo dwarfing or malforma- 
tion). The neutralizing titer of each serum against homologous or heterologous virus was determined by the last 
serum dilution that protected 50% of the embryo and calculated by Reed-Muench method”. In addition, the 
cross neutralization r-values between strains were calculated by the method described by Archetti and Horsfall°". 
Antigenic (serotype) difference between two given strains was defined as follows, when r: 70-100%, same sero- 
type; r: 33-70%, different subtype (minor); r: 11-32%, different subtype (major); r; 0-10%, different serotype’. 


IBV infections in SPF chickens. For survival evaluation after IBV infection, 1-day-old SPF chicks were 
randomly divided into three groups and intranasally challenged with 10° EIDsy of IBV 3575/08, 2575 WT and 
2575 AT, respectively (n=6 per group). Another six chickens received PBS as naive controls. For the clinico- 
pathological assessment, 2-week-old SPF chickens were used. Following intranasal challenge with 10° EID.) of 
IBV 3575/08, 2575 WT and 2575 AT, chickens were monitored on a daily basis. The clinical scores of IBV were 
interpreted according to the methods described by Avellaneda et al.*°. The clinical signs were evaluated as: 0 = no 
clinical signs; 1 = slight lacrimation, slight shaking, watering feces or tracheal rales; 2=lacrimation, presence of 
nasal exudate, depression, water feces, apparent sneezing or cough; 3 = same as 2 but stronger with severe watery 
feces; 4= death. Throat and cloacal swabs were collected at 1, 2, 3, 7, 14, 21 dpi for virus detection. In addition, 
lungs and kidneys were also harvested from euthanized chickens at 3, 7, 21 dpi for virus detection and histopatho- 
logical examination. Sera were collected at 14 and 21 dpi to test the IBV antibody titer. For PCR array analysis, 
three chickens at 24 and 48 hpi from each challenged group were euthanized to harvest lungs. Four chickens 
received PBS as naive controls. For all chicken experiments, each virus-challenged group of chickens was reared 
in separated rooms, and utensils were separately handled to avoid any cross contamination. 


Viral detection by RT-PCR. The viral nucleic acid extraction kit (Geneaid) was used to extract the viral RNA 
from swabs samples and tissue homogenates. One step RT-PCR was performed in a total reaction volume of 25 il 
containing 0.5 tl Phire hot start II DNA polymerase (Thermo), 5 11 5X buffer, M-MLV RT (200 U/l) (Invitrogen), 
1p1 RNaseOUT (40 U/l) (Invitrogen), 0.5 11 501M forward primer (5’-GGTAGYGGYGTTCCTGATAA-3’) and 
reverse primer (5’-TCATCTTGTCRTCACCAAAA-3’) (Tri-I Biotech) for a 618 bp fragment of IBV N gene 
detection, 3 1L 2.5mM dNTPs (GeneTeks BioScience), 2.5 ul RNA, and 10.5 il nuclease-free water. RT-PCR was 
performed at 40°C for 30 min, 98°C for 3 min, followed by 35 cycles of denaturation at 98°C for 30s, annealing 
at 58°C for 30s, and extension at 72°C for 30s, and the final extension was performed at 72°C for 3min. PCR 
products were then analyzed on 2% agarose gels. 


ELISA. A flat-bottomed microplate (Nunc, Denmark) was coated overnight with purified IBV antigen (100 ng/well) 
with coating buffer (15 mM Na,CO, and 35 mM NaHCO; , pH 9.6). The wells were washed three times with PBST 
(0.1% (v/v) Tween 80 in PBS) and blocked with 5% (w/v) skim milk (BD Difco, Sparks, MD) at room temperature 
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for 1h. After blocking, serial diluted chicken serum was added to each well (100 l/well) and incubated at room 
temperature for 1 h. Following three washes, 100 il of 1:2,000 diluted peroxidase-conjugated affinipure rabbit 
anti-chicken IgG (H+ L) (Jackson ImmunoResearch, West Grove, PA) was dispensed into each well and incu- 
bated at room temperature for 1h. After three additional washes, the wells were then incubated with 100 1l of 
SureBlue Reserve TMB Microwell Peroxidase Substrate (KPL, Gaithersburg, MD) for 10 min in the dark. The 
reaction was stopped by the addition of 100,11 of TMB stop solution (KPL). The optical density (OD) at 450nm 
was read using an automated plate reader (Thermo). 


Histopathological examination. Chicken tracheas, lungs, and kidneys were collected and fixed in 10% 
formalin (Sigma, MO) at room temperature for at least 48h. Fixed tissues were embedded in paraffin wax, and 
cut into 61m sections. Sections were then routinely processed and stained with hematoxylin and eosin (H&E). 
Microscopic lesions resulted from IBV infection were examined by light microscopy. 


PCR array. To analyze the innate and adaptive immune responses in chicken after infection, RT? Profiler 
Chicken Innate & Adaptive Immune Responses PCR Array (PAGG-052ZG; SABiosciences) was used following 
the manufacturer's protocol. Total RNA in chicken lung was extracted using TriSolution Reagent Plus (GMbiolab, 
Taipei, Taiwan) following the recommended protocol for tissue samples and quantitated using spectrometry. To 
avoid the contaminations of genomic DNA, 1 jg total RNA was treated with 1 U DNase I (Thermo) for 30 min at 
37°C, then 1 ul EDTA was added and incubated 10 min at 60°C to inactivate the DNase. cDNA was synthesized 
using the RT? First Strand kit (SABiosciences), mixed with 2x RT? qPCR Master Mix (SABiosciences), and then 
loaded onto the 384-well PCR array plate for amplification on the LightCycler 480 PCR system (Roche). The steps 
of the cycling program were 95°C for 10 min followed by 40 cycles of 95°C for 15s and 60°C for 1 min. Data were 
analyzed using the online software provided by SABiosciences (http://pcrdataanalysis.sabiosciences.com/pcr/ 
arrayanalysis.php). 


Real-time RT-PCR. Total RNA in chicken lung was extracted as described above. Real-time RT-PCR was 
performed with QuantiNova Reverse Transcription kit (Qiagen, Alameda, CA) and QuantiNova SYBR Green 
PCR kit (Qiagen) using previously described primers that target the chicken TLR3, TLR4, and 28S™. All reactions 
were set up in duplicate on the iQ5 real-time PCR system (Bio-rad). The steps of the cycling steps were 95°C for 
2 min followed by 40 cycles of 95°C for 5s and 60°C for 10s. Obtained cycle threshold (Ct) values were normal- 
ized to 28S. The relative expression of TLR3 and TLR4 (fold change of naive control) was determined using the 
2-44¢ method®. 


Statistical analysis. Data were analyzed by unpaired t tests or ANOVA followed by Dunnett’s multiple com- 
parison tests using GraphPad Prism (GraphPad Software, San Diego, CA). The p values smaller than 0.05 were 
considered significant. 
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